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ABSTRACT 

The journey of the Sun through space carries the solar system through a dynamic 
interstellar environment that is presently characterized by a Mach ~ 1 motion be¬ 
tween the heliosphere and the surrounding warm partially ionized interstellar cloud. 
The interaction between the heliosphere and interstellar medium is an evolving pro¬ 
cess due to variable solar wind properties and the turbulent nature of the interstellar 
cloud that surrounds the heliosphere. Frisch et al. presented a meta-analysis of the 
historical data on the interstellar wind flowing through the heliosphere and concluded 
that temporal changes in the ecliptic longitude of the flow direction with time were 
statistically indicated by the data available in the refereed literature at the time of 
that writing. Lallement & Bertaux disagree with this result, and suggested, for in¬ 
stance, that a key instrumental response function of IBEX-Lo was incorrect and that 
the STEREO pickup ion data are unsuitable for diagnosing the flow of interstellar neu¬ 
trals through the heliosphere. In this paper we first show that temporal variations in 
the interstellar wind through the heliosphere are consistent with our knowledge of the 
very local interstellar medium. The statistical analysis of the helium wind data is re¬ 
visited, and a recent correction of a typographical error in the literature is incorporated 
into the new fits. With this correction, and including no newer IBEX results, these 
combined data still indicate that a change in the longitude of the interstellar neutral 
wind of A = 5.6° ± 2.4° over the past forty years remains statistically likely, but an 
constant flow longitude is now statistically possible. Other scenarios for the selection 
of subsets of these data used in the fitting process produce similar conclusions. We 
show that the speculations made by Lallement & Bertin about the IBEX instrumental 
response function are incorrect, and that their other objections to the data used in the 
meta-analysis are either incorrect or unproven. Further investigations of the historical 
interstellar wind data and continuing analysis of additional IBEX data may provide a 
more definitive result on the stability of the flow of interstellar gas through the helio¬ 
sphere. 

Subject headings: ISM: magnetic fields, heliosphere— solar system: general 


1. Introduction 

The solar system is traveling through a dynamically evolving interstellar medium that leads 
to variations in the interstellar wind through the solar system over geologically short time-scales 
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( Frisch et~aD 2011: Fri sch & Mueller il2013 ). The Loc al Interstellar Cloud (LI C) now around the he¬ 
liosphere is a low density and partially ionized cloud (Slavin & Frisch 2008). The relative motions 
of the Sun and LIC creates a wind of interstellar gas and dust through the heliosphere whose veloc¬ 
ity was first defined by measurements of the fluorescence of solar 584A photo ns off of interstellar 
He° in the gravitational focusing cone downwind of the Sun (Welle r & Meier 197 4). The absence 
of interstellar H° at the LIC velocity toward the two nearby stars a Cen (1.3 pc) and 36 Oph (6.0 
pc) indicate t hat the Sun will emerge from the LIC in t he upwincQ direction anytime within the next 
3,500 years ( Lallement et al. 1995 : Woo d et al. 2000 ). and that it entered the LIC during the past 
10,000 years (Frisch ll994l) . Variations in the interstellar wind direction over historical timescales, 
of about 40 years for the space age, suggest that interstellar neutrals, and the pickup ions that form 
from them after ionization, provide the best diagnostics of the longitude of the interstellar wind 
flowing through the heliosphere. The dynamic characteristics of the surrounding LIC suggest that 
short-term variations in the interstellar wind direction may exist. 

In the paper titled ’’Decades-Long Changes of the Interstellar Wind Through Our Solar Sys¬ 
tem”, Frisch et al. (2013, hereafter Paper I) have presented a meta-analysis of historical obser¬ 
vations of the interstellar He° wind data and found that the data indicated a statistically likely 
shift in the wind longitude of ~ 6.8° over the past forty years, and that a constant flow lon¬ 
gitude is statistically unlikely. The historical data set included in situ measuremen ts of the he¬ 
lium wind by IBEX ( McComas et al. 2012 : Mobius et al. 2012 : Bzowski et al. 2012 ) and Ulysses 
(I Wittel 120041) . measurements of pickup ions in the downwind gravitational focusing cone by the 


ACE-SWICS, STEREO-PLA STIC, and MESSENGER-FIPS (IGloeckler et al.ll2004l: Drews et al. 


20121: iGershman et al.1120131) . the upwind crescent by STEREO-PLASTIC (IDrews et al.l 120121) . 


interstellar He° 

in the inner heliosphere (iWeller & Meier 

1974 

Aiello 

1978 

; l 

Velio et al. 1979; 

Weller & Meier 

1981; 

Dalaudier et al. 1984; 

Vallerga et al. 

2004; 

Nakagawa et a 

.2008 

). This was 


the first meta-analysis that searched for possible historical variations in the direction of the wind 
of heavy interstellar neutral atoms through the heliosphere. 

The original study in Paper I was restricted to heavy atoms and did not include longitude infor¬ 
mation from the velocity vectors of either interstellar H° or interstellar dust grains flowing through 
the heliosphere. Both sets of particles are strongly filtered on a time-variable basis, beca use of the 
solar magnetic activity cycle, throughout the heliosphere and heliosheath regions (IRinken & Fahi 


1983; 

Rucinski & Bzowsk 

1 1996; 

Sau 

let al. 2013; 

2012; 

Sterken et al. 

2012, 

Strub et al. 

2014b- 


et al.ll20 13 : Grogan et~alll996 : Landgraflhoool: Slavin et al. 


'The terms “upwind” and “downwind” are referred to the direction of the heliosphere nose, which in turn is defined 
by the direction of the interstellar He° wind that flows into the inner heliosphere. 
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In a paper_”On the Decades-Long Stability of the Interstellar Wind through the Solar Sys¬ 


tem”, Lai lenient & B e rtaux (120141 hereafter LB 14) disagree with the conclusions in Paper I. The 
different viewpoints presented in Paper I and LB 14 fall into two categories, the most interesting 
of which is the underlying science question that is encapsulated by the differences of the two titles 
”On the ’’Decades-Long Changes of the Interstellar Wind Through Our Solar System” compared 
to ’’Decades-Long Stability of the Interstellar Wind through the Solar System”. The second area of 
disagreement centers on the selection of data used in the analysis, and the IBEX instrumental pa¬ 
rameters. LB 14 favor the interstellar He° flow direction determined in Mob ius et ah (l20041) . where 
measurements of the interstellar He° wind found by Ulysses, ACE, EUVE, and Prognoz 6 were 
combined to find the best-fitting direction for the interstellar wind. The groundwork for under¬ 


standing the interstellar helium wind is presented in the set of papers that accompany Mobius et ah 
J2004i) Pl 


In this paper we first revisit the conclusions of Paper I, with the proviso that new data on, and 
better modeling of, the interstellar wind are creating a rapidly changing scientific data base upon 
which our conclusions rely. Data showing that the interstellar wind could change over historical 
times are discussed, and these data are consistent with inhomogeneities in the LIC over spatial 
scales comparable to the heliosphere dimensions Lj2]). New statistical analyses are performed on 
the historical data (IQ, utilizing a new published correction for the uncertainties on the NOZOMI 
interstellar He° wind direction ( 114.3.31) . A temporal change in the wind longitude is found to be 
statistically likely, but a constant wind direction also becomes statistically acceptable although less 
likely than a variable direction. LB 14 objected to the use of the IBEX-Lo data in recovering a 
longitude for the He° flow that differed from the value found by Ulysses. Methods for retrieving 
the flow longitude from IBEX-Lo data are discussed in 1 14.11 with further details on the instru¬ 
mental response function and properties in 1 14.1.31 and Appendix |Bj Cloud temperature plays an 
important role in the uncertainties of the parameter space constrained by the IBEX-Lo data; the 
interstellar picture of the LIC temperature is discussed in Appendix lAl The LB 14 comments on the 
STEREO data are discussed in 114.21 where new data are mentioned that contradict the hypothesis 
of significant azimuthal transport of pickup ions inside of 1 AU, and in Appendix O 


2 The wind direction in this study needs to be corrected by 0.7° to account for an unrecognized difference in 
coordinateepochs in the contributing data that led to an error of about 0.7° in the wind direction ( Frisch et al. 20091: 
Lallement et al.lEoiO ). 
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2. Possible spatial inhomogeneities in LIC 
2.1. Turbulence and edge effects 


The Sun moves through the surrounding interstellar cloud at a relative velocity of ~ 5 AU 
per year, so that the 40-year historical record of the interstellar wind velocity sampled interstellar 
scale lengths of ~ 200 AU. There is no ad hoc reason that the LIC should be homogeneous and 
isotropic over such small spatial scales. The mean free path of a thermal population of LIC atoms 
is ~ 330 AU@ which is larger than the distance traveled through the LIC during the past 40 years. 
Collisional coupling between atoms will break down over scales smaller than the mean free path, 
allowing the formation of eddies that perturb the gas velocity over heliosphere scale-lengths. Paper 
I suggested that a shift in the He° flow direction could be due to the nonthermal turbulence of eddy 
motions. 


A second source of small -scale turbulence in the LIC is Alfvenic tur bulence. The LIC inter¬ 
stellar magnetic field (ISMF) (McComas et al. 20091: Funsten et al. 20131) affects clumping of the 
charged particles over scales smaller than the collisional mean free path. The LIC magnetic field 
strength is ~ 3 jiG, based on the plasma pressure of energetic neutral atoms (ENAs) in the upwind 
direction and the magnetic distortion of the heliotail (ISchwadron et aU20111). and the eq uipartition 
of energy between magnetic and thermal energies in the LIC (|Slav in & Fri s ch 2008) . For field 
strengths of 3 /iG and larger, the Lorentz force will couple protons tightly to the ISMF, over ~ 375 
km scales for a gas temperature of 7000 K. This tight coupling acts as a conduit for magnetic turbu¬ 
lence into the LIC that affects the neutral population observed by IBEX through charge exchange, 
and would perturb the neutral particle velocities over spatial scales smaller than the collisional 
scale. 


The location of the LIC in a de celera ting group of clouds that appears to be a fragment of 
the Loop I superbubble shell (Frisch et al. 201 F) wo uld promote turbulence in the local ISMF. If 
the surrounding decelerating flow of ISM (Frisch et al. I2002L 20111) consists of clouds in contact 
with each other, then macro-turbulence may form in the LIC at the acoustic velocity of about 8 
km s' 1 . The angle between the LIC ISML and the flow of LIC gas through the local standard 
rest (LSR) is 86° =b 14°, indicating that the velocity and magnetic field vectors are perpendicular 
(Lri sch & Schwadron 2014). Depending on the relative thermal and magnetic pressures, compres¬ 
sive turbulence may form. Magnetic turbulence may propagate along the bridge of po larized dust 
that extends from near the Sun out to the North Polar Spur region (Lrisch et al. 2014). Magnetic 
disturbances will propagate along the ISML at the Alfven speed of ~ 25 km s' 1 (using an electron 


3 The mean fre e pa th in the LIC is dominated by charge-exchange and induced dipole scattering between H° and 
H + atoms ( Spangler et all20U ). 






































- 6 - 


density of 0.07 cm 3 , Slavin & Frisch 20081) . organizing additional disruption of particle veloci¬ 


ties over scale sizes smaller than the mean free path. 


It is important to note that the thermal distribution of the He° atoms smooths the time-interval 
sampled by the population observed at 1 AU because of the different propagation times of fast and 
slow atoms through the heliosphere. A Maxwellian gas at 7000 K has a Doppler width of ~ 5 
km s _1 , which is a significant percentage of the bulk LIC motion, causing faster particles to catch 
up to slower particles in the incoming He° flow. The Maxwellian distribution will not apply over 
microscales less than the 330 AU mean free path of charge-exchange, and in principle the mean 
free paths of the higher energy atoms will be larger than 330 AU anyway. 


Upper limits on an interstellar Lyct component at the LIC veloci ty toward 36 Oph, in the 
upwind direction, show that the Sun is near the LIC edge (Wood et al. 2000). The high level of 
LIC ionization that is inferred from photoionization models requires that ~ 22%, 39%, 80%, and 
19 % of th e h ydrogen, helium, neon, and oxygen atoms at the heliosphere boundaries are ionized 
(Slavi n & Fri sch 120081. SF08). A thin conductive interface between the warm LIC and hot plasma 
of the Local Bubble maintains these ionization levels (SF08). A sharp density decrease and velocity 
increase is expected at the edge between the LIC and conductive interface. The unknown distance 
to the upwind edge of the LIC allows that such an interface is close to the heliosphere. However, 
no direct observational evidence of this interface is yet found. 


2.2. Supersonic collisions between clouds 


Disturbances due to the propagation of shocks through the cloud are also possible sources of 
kinematical structure in the LIC. The LIC belongs to a decelerating flow of ISM through space 
that may create dynamical interfaces in the interaction regions between clouds in physical contact 
(Frisch et al. 2 0111) . If the LIC is in physical contact with either the Blue Cloud (BC) in the 
downwind direction, or the G-Cloud (GC) in the upwind direction, 0 the shocks that form at the 
supersonic interfaces will drive compressional disturbances into the LIC that perturb the density 
and velocity of the gas particles. 

A supersonic interface would form between the LIC and GC if they are in contact. The angle 
between the LIC and GC LSR velocity vectors is 17.2°, and the relati ve ve locity is 14.6 km 1 . 
The G-cloud has a temperature of 5500 ± 400 K (IRedfield & Linskvl 20081) . corresponding to an 


4 An independen t sor ting of local interstellar velocity absorption components into clouds is provided by 


Redfield & Linskv (2008) and we use that cloud-naming scheme here. 
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acoustic velocity of ~ 7.7 km s \ indicating a shock strength of Mach~ 1.9 between the cloudsjfl 


A supersonic interface would also form between the LIC and the BC if they are in contact. 
These clouds have LSR velocities of 15.3 ± 2.5 km s” 1 and 7.1 ± 2.3 km s -1 respectively, and an 
angle between the two vectors of 49° (cloud LSR velocities are from Frisch & Schwadron 2014). 
The LIC sound velocity is 8.6 km s’ 1 , for a 7000 K perfect gas with mean molecular weight 
1.29. The ~ 11.9 km s’ 1 relative velocity indicates a supersonic interface of Mach number ~ 1.4 
between the two clouds if they a re in direct contact. The LIC is warmer than the BC, which has a 
temperature 3000 + 2 oqo K (iHebrard et al.lll999h . A simple temperature gradient between the three 
clouds does not exist. In the presence of the interstellar magnetic field, the shock transitions will 
be modified by the unknown angle between the magnetic field and the normal to the cloud surface. 


3. Refitting the historical data on the interstellar wind longitude 


A new statistical interpretation of the historical data (Table 1 in the online supplementary 
material, OSM, o f Pape r I) is required by a typographical error in the original NOZOMI publication 
( Nakagawa et alJ 2008). We first perform new fits to the same set of data used in the Paper I 
using th e corre cted uncertainty for the He° flow longitude determined from the NOZOMI data 
( Nakagawa et al.ll2014l 1)4.3.31) . and next for an alternate selection of the fitted data. 


The results of the refitting are presented in Table [2] The first result, heading A, shows the 
original fit from Paper I. The updated fit using the corrected NOZOMI uncertainties is listed under 
heading B, fits 3 and 4, of Tabled and plotted in Figured] The slope of the revised fit is 0.14±0.06, 
which implies a directional change of 5 A = 5.6° ± 2.4° over the past 40 years (instead of the 
6.8° ± 2.4° in Paper I). The reduced y 2 of the fit is 0.7, and the p-value of 0.82 is larger than the 
0.05 value required to pass the p-test [j so that this result is highly likely. The small value for y 2 
suggests that the uncertainties may be overestimated by about ~ 20%. The corresponding fit for 
the assumption of a constant flow longitude gives A = 75.0 ± 0.3 (fit 4), which has a p-value of 


5 Using instead the heliocen tric v elo cities for th e LIC and GC in the upwind direction, and now-obsolete physical 
parameters for both clouds, Grzedzielski & Lallement (1996) have predicted that a supersonic interface is present 
between the LIC and the G^cloud, with the LIC corresponding to the post-shock gas because it is warmer than the 
G-cloud. Frisch _& York ( 199 li) also modeled the local shock structure. 

6 As pointed out in S5 of the OSM in Paper I, the probability of taking a result y 2 more extreme than the observed 
value of Xobs * s given by the p-value that equals the minimum between the two probabilities, P(0 < y 2 < y 2 bs ) 
and P(y 2 bs < y 2 < oo). A null hypothesis associated with a p-value smaller than the significance level of 0.05 
is typi c ally rejected. Ad d itional informa ti on about t he sta tistical techniques used in the fits can be found 
( 2014 ). Livadiotis ( 2007 ). Deming ( 1964 ). Lvbanon h 1984 ) 


Livadiotis 
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0.014 and is now a statistically acceptable fit, but with a lower likehood than a fit with a variable 
flow direction. This new fit to the data using the corrected NOZOMI uncertainties should replace 
the fit that is presented in Paper I. 

In order to test the robustness of the fit to these data, we have performed a second fit with four 
additional modifications made to the input data set: 


1. The uncertainties assigned to the STP 72-1, Mariner 10, and SOLRAD1 IB flow directions in 
Paper I represented the combined longitude and latitude uncertainties reported in the original 
publications ( E14.3I) . The new fit uses only the longitude uncertainty. It is assumed that the 
uncertainties for longitude and latitude are normally distributed, so that the uncertainty used 
in Table 1 of Paper 1 can be reduced by a factor of 1.414. 


2. The IBEX and Ulysses data are treated equivalently in the sense that the two independent 
data sets corresponding to the 2009 and 2010 seasons are included separately,_jather than 
combined as in Paper I. The longitudes for the individual seasons from Bzowski et al. (2012) 
are used. 


3. The uncertainties for the IBEX longitudes are based on the IBEX data alone, without re¬ 
course to astronomical data on the LIC. The IBEX temperature uncertainty range of 4400- 


8200 K ( (14.1.2b can be used to set longitude uncertainties of ±4.2° (Fig. 1 of McComas et al. 


2012). In Paper I a LIC temperature based on astronomical data was used, however the astro¬ 


nomical data reveal large variations in the LIC temperature toward different stars (Appendix 

m. 


4. The Prognoz 6 data are omitted from this fit because both the geometric fit and the fit ob¬ 
tained from modeling the data are no longer considered valid by coauthors of the original 
papers (Section 14.3 .2b . 


Three fits were performed on this adjusted data set, a linear fit, a fit assuming an constant flow 
longitude, and a parabolic fit (heading C, fits 5, 6 and 7). The goal of the fits is to test whether a 
temporal variation in the flow direction over time is statistically more likely than a constant flow 
direction over time. The linear fits listed in Table |2] are highly likely and always more likely than 
a fit with a constant longitude over time. The low x 2 value of the parabolic fit suggests that the 
limited set of available data does not justify a second-order fit. Note that the finding that a linear fit 
to the data is statistically more likely than a constant or parabolic fit is not equivalent to a statement 
that a change in longitude must have occurred in a linear fashion. 

The linear fit is statistically the most likely fit and gives a temporal variation of 5.6° ± 2.4° 
over the past 40 years (the fit is shown in Figure [0 right). This result is close to the result in Paper 
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I. However while a constant longitude could be rejected in Paper I, the increase in the NOZOMI 
uncertainties makes a constant longitude over time acceptable but less likely than a variation of the 
longitude over time[Livadiotis (2014). 


The parabolic fit is included to show that a fit with a higher order polynomial yields a result 
that is less satisfactory than a fit with a lower order polynomial. As the order of the polynomial in 
the fit increases, the x 2 should in principle decrease until it is < 1 or it is x 2 ~ 1- The parabolic fit 
under heading C gives 0.6, which tells us that the best statistical description of the temporal 
variation in the flow direction is given by the linear fit. 


New reductions and modeling of the U l ysses data have appeare d since the publication of Pa 


per I ( Bzowski et all2014 : Wood et ap2014 : Katushkina et al. 


2014). These results find a He° flow 


longitu de for the Ulysses data collected during the 2007 polar pa ss that is close to the v alue of lwittel 
( 2004 ). but with a somewhat higher temperature (7500^2ooo K, Bzowski et al. 2014), The differ 


ences between the He° flow vector found from the IBEX 2009-2010 seasons of data (M obius et al. 


20121: iBzowski et al.ll2012c iMcComas et al.ll2012t) . and the Ulysses vector, are not yet explained. 


For t he latest seasons of IBEX observations (2012-2014, IMcComas et al.1 1 201 4c iLeonard et al. 


2Q14), the results from spacecraft pointings in the ecliptic suggest a di fferen t portion of the cou¬ 


pled 4-dimensional “tube” of possible solutions presented by McComas et al. (2012j) with similar 
He° velocity vectors to Ulysses but requiring a significantly warmer cloud; other new IBEX point¬ 
ing directions, at different angles to the ecliptic plane, are not consistent with the Ulysses data for 
reasons not fully understood. Comparisons of Ulysses and IBEX results are a study-in-progress 
and we do not include results presented after 2013 in the refitting of the data in this paper. 


4. Data used for detecting temporal variations in the neutral helium wind 


LB 14 have argued that most of the data incorporated into the meta-analysis of Paper I were 
either incorrect or were used incorrectly. These data, all of which were taken from the refer¬ 
eed scientific literature, fall into three groups: remote measurements of the backscatter ed so- 


lar 584A emission from interstellar He° in the inner heliosphere ( 

Weller & Meiei 

1974 

Meier 

1977; 

Broadfoot & Kumar 

1978; 

Aie 

lo 1978j; Aiello et al 

1979; 

Weller & Meier 

1979L 

1981; 

Dalaudier et al. 1984; Flvnn et al 

1998 

; Vallerga et al. 2004; 

Lallement et al. 2004; Nakagawa et al. 

2008, 

2014); in situ measurements of the pickup-up ions that are sampled by spacecraft travers- 


ing the gravitational focusing cone and upwind crescent (Gloeckler et al. 2004; 

Drews et al. 2012: 

Gers 

iman et al. 2013), and in situ measurements of interstellar He° atoms (Witte 

2004; Mobius et al. 

2012 

Bzowski et al. 2012; McComas et al. 2012; Bzowski et al. 2014). 



The main LB 14 arguments are summarized in Table Q] along with our comments on those 
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arguments. Most of the arguments are speculative. The LB 14 hypothesis about an underestimated 
IBEX-Lo instrument response function is wrong ( ^4.1.31 ! fBl) . The better-supported argument that 
the STEREO pickup ion data only provide upper limits because of the azimuthal transport of 
pickup ions inside of 1 AU has since been shown to be i ncorre ct by the 0.1 AU parallel mean-free- 
path of newborn pickup ions found bv lGershman e t al. (2014) from MESSENGER data ( 44.21) . 


4.1. Evaluation of He° flow longitude from IBEX measurements 

4.1.1. Functional dependencies ofHe°flow direction on peak neutral fluxes 


The trajectories of neutral interstellar helium atoms through the heliosphere depend on the 
four parameters describe the interstellar He° gas at “infinity” beyond the influence of the helio¬ 
sphere, the gas flow longitude, Aism,oo> latitude /3 ism,oo, speed V IS m,<x» and temperature T ISM]0O . 
In the parameter space measure d by IBEX-Lo, these parameters are in terde pendent with c ouple d 
uncertainties. As discussed by iMobius et al.l (120121) . iLee et al.1 (120121) and iMobius et al.l (120141) . 


since IBEX is nearly a Sun-pointed spinner, the determination of the interstellar neutrals (ISN) 
inflow vector is tied to the IBEX-Lo observation of the flux maximum that is in a direction that 
is perpendicular to the Earth-Sun line. Then, as shown in Figure O the observed flux maximum 
is directly related to the bulk flow vector at 1 AU, in the plane of the ISN bulk flow trajectory. 
Here Vism,oo and Ai S m,oo are uniquely connected by the hyperbolic trajectory equation (eqn. Q]) 
and the relation between the true anomaly (the angle between the observed velocity and the 
flow vector at infinity) of the trajectory, the observer location when obtaining the peak A pea k, and 
inflow direction Ai SM ,oo, so that 6> oo =Ai S M,oo+180°-A peak is: 


VlSM.oo — 


\GM S 

( _1 ill 

r E 



1/2 


( 1 ) 


Without further information, this constitutes a functional relationship of Vism,oo(Aism,oo)- The 
location of the flow maximum observed by IBEX along the Earth orbit A pea k was determined to 
be 130.6° ± 0.7° in (Mobius et al. 2012 ). which constrains the functional relationship in eqn. 1 to 
within < ±0.7° in Aism,oo and < ±0.5 km s -1 in V IS m ! 00 , although a large range of values along 
the function satisfy the observed location of the interstellar gas flow maximum at 1 AU nearly as 
well. 

Simultaneously, the ISN flow direction in latitu de .^Ts vT.-y.CArsM.or,) and the ISN temperature 
Tism,oo(Aism,oo) are determined in the IMobius et al.l (2012) analysis from the ISN flow peak lo¬ 
cation in latitude and its angular width in latitude, o the ob serva tion of whi ch is shown in Fig¬ 
ure [3] again varying with the ISN flow longitude (also see Mobius et al. 2014). As indicated in 
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the left panel of Figure [2 the observed angular distribution of the incoming interstellar neutrals 
constitutes a local Mach cone of width vti/Vism,iau> based on the bulk flow Vism,iau and the 
thermal velocity i; Th of the distribution at 1 AU. Therefore, the observed angular width fixes the 
ratio T^MjCc^/Vkmpai/* and the derived temperature depends on the derived ISN bulk flow speed 
V ism.oo? the only variable in Vism,iau> and thus on Aism,oo- Because the statistical uncertainties of 
the observed peak location and width of the ISN flow are only 2% and 1%, respectively, the com¬ 
bined uncertainties of the two relations are mainly determined by the uncertainties in relation (1) 
and potential systematic effects. We will address such effects below, in particular, on the resulting 
temperature. 


As shown in Figure |4l the latitude V’Peak of the flow at 1 AU in the solar rest frame is connected 
to the ISN flow latitude /3 ism,oo at infinity through the true anomaly of the bulk flow as 


tan( , 0. 


Peak ) 


tan (/3 ism ,oo) 


tan(/3isM ,C 


sin(6> 00 ) | | sin(AisM,oo + 180° - A Peak ) | 


( 2 ) 


and thus again to the ISN inflow longitude Aism.oo- These relations constitute a degeneracy in the 
determination of V ISMi00 , Asm.oo, and T ISMi00 as a function of Aism.oo- However, relation (2) also 
provides a tool to break the degeneracy using IBEX latitude distributions combined for different 
positions of IBEX as the observer varies in longitude A obs- Then, A obs replaces A pea k in equation 
(2) as a variable and the resulting relation ^(Aism.oo) depends distinctly on Ai S m,oo- The comparison 
of ISN flow distributions observed over entire 2009-2010 observatio n seas ons with the predicted 
relation (2) as a function of Aism,oo and respective x 2 -minimization in Mobius et al. (2012 ) resulted 
in the reported optimum inflow longitude Ajsm ^=79.0° (+3.0°, —3.5°) with a similar result in 


Bzowski et al. (2012). iBzowski et al.1 (120121) used a x 2 analysis to find an uncertainty r ange for 
the interstellar flow longitude of ~ ±4° (blue contour in Fig. 22 of IBzowski e t al . 2012). These 
uncertainties were the basis for the bounding values of the resulting parameter range from the 
IBEX observations used in McComas et al. (2012). 


In contrast, Lallement & Bertaux (2014, in top box of Fig. 2) maintain that the variation of 
the observed fluxes as a function of observer longitude A obs, with a substantially wider maximum 
than observed with Ulysses GAS (Witte 2004), may have led to the different inflow longitude, 
which may be affected by instrumental influences. To support their claim, they compare the ob¬ 
served IBEX count rate distributions with simulated distributions for the Witte (2004) parameter 
set produced by Bzowski et al. (2012) from the best-fitting x 2 analysis. However, the difference 
in the distributions between these two sets is totally dominated by the larger angular widths of the 
count rate sequences measured by IBEX, which leads to the larger derived temperature, and not 
by unaccounted-for data losses. As shown by Lee et al. (2012), the widths of the distribution in 
latitude and longitude are tightly coupled and represent the ISN temperature. 

To illustrate the basis for the inflow direction in the IBEX analysis, Figure [5] shows the peak 
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location and width of the ISN flow distributions in latitude as simulated for the observation periods 



found by Witte (2004), which stands out as much narrower than that from IBEX. However, the 
visible differences between the two parameter sets along the IBEX tube are, first, a steeper variation 
of the peak location as a function of orbit for Aism,oo as favored by Bzowski et al. (2012) than for 
the value of Witte (2004) and, second, a somewhat stronger increase in width towards the_early 
orbits. To arrive at a value for Aism,oo, solely the first effect has been used by Mobius et al. (2012). 
In relation to Aism,oo the ^'-minimization in the comparison between observed and simulated flow 
distributions used by Bzowski et al. (2012, see Fig. 22) is sensitive to both effects as illustrated in 
Figure [5] 


This leaves only the argument of LB 14 that the width of the angular distribution, which now 
solely influences the derived temperature, could be brought into agreement with the Ulysses value 
by considering the influence of a strong and unaccounted-for instrumental effect in the IBEX data 
reduction. As derived from extensive tests with the actual IBEX data system hardware simulator 
and Monte Carlo simulations of realistic data throughput of a wide range of particle rates, the actual 
effect of data suppression due to limited transmission bandwidth on the width of the observed 
distributions is < 2% and thus is negligible (see Appendix |B] including Figure [6]). 


LB 14 make the point that the degeneracy between the longitude and velocity of the He° distri¬ 
bution function allows the longitude of the peak count rates of IBEX data to also correspond to the 
expectations for the Ulysses flow vector, providing there is an additional suppression of count rates 
(detector “dead-time”) that was neglected in the IBEX data flow)] In support of their hypothesis 


that the helium vector of 


Witte (2004) is correct, LB 14 took the orbit-by-orbit IBEX data, and the 


simulations where Bzowski et al. 


(120121 Fig. 18) compared the IBEX data with the Ulysses flow 
parameters. Those rates and data simulatio ns were rep lotted for two cases (Fig. 2 of LB 14). One 
line on the plot showed the results of the B z owski et a]J (2012) analysis, while the other line plotted 
an adjusted count-rate based on the LB 14 hypothesized 7 ms dead-time correction for the IBEX- 
LO instrument. They argued that the Witte (2004) flow direction agreed with the IBEX count rates 
for their hypothesized dead-time correction^ This argument to reject the IBEX analysis because 


7 More specifically, they state that “We have examined the consequences of dead-time counting effects, and con¬ 
clude that their inclusion at a realistic level is sufficient to reconcile the [IBEX 2009-2010] data with the old [Witte 
Ulysses] parameters, calling for further investigations.” 

8 LB14 state: ”...it is premature to use the new longitude and velocity [of IBEX] before the question of the dead-time 
correction is fully settled.” 






















-13- 


of a hypothesized 7 ms dead-time correction is not valid. 


The effect that influenced the IBEX data flow is a known finite data processing time, not a 
sensor dead-time, and the types of events mostly responsible for this effect were unanticipated elec¬ 
trons, which do not alter the shape of the observed interstellar flow distribution significantly. Also, 
these are excluded by a new instrument mode that started mid-2012. This operational modification 
has eliminated the count rate reduction reduction effect almost entirely. Comparison between the 
old and new operational modes has demonstrated that the finite data processing time does not sig¬ 


nifican t ly affect the count rate distri butions used to reach th e conclusions reported i mMobius et al. 


(120121) . iBzowski et al.1 (120121) . and McComas et al.1 (120121 Appendix |B]). In particular, the data 
processing time did not influence the deduced He temperature, which LB 14 tried to reconcile by 
exploring this effect to an extent that was already excluded by Mobius et al. (2012). The instru¬ 
ment response function of IBEX-LO does not contribute in any significant way to the uncertainties 
of the interstellar He° flow direction derived from the IBEX data. This is discussed in more detail 
in Appendix |B] 

In summary, the inflow longitude parameter hinges mostly on the evolution of the peak loca¬ 
tion in latitude as a function of observer longitude in the IBEX observations. While the temperature 
in Mobius et al. (2012) is obtained by fitting the height and width of the count rate distributions, 


these observables were ignored in LB 14. The temperature derived from the IBEX observations, 
for the case where one assumes the inflow velocity of the Ulysses results (Witte, 2004), is 30% 
higher than the temperature reported from Ulysses GAS observations of He°. Systematic instru¬ 
mental effects larger than 2% in the IBEX observations, such as an erroneous instrument response 
function, can be ruled out as the reason for the difference between the Ulysses and IBEX He° flow 
directions ( also see Appendix®. 


4.1.2. IBEX results for 2009-2010 seasons of data 
IBEX-Lo directly samples the trajectories of interstellar He atoms at 1 AU after they have 


reach the 1 AU orbit of IBEX ( 

McComas et al. 

20oi 

Aisclier et al. 

2009b 

Mobius et al. 

2009; 

Fuselier et al. 

2009a; 

Mobius et al. 

2012; 

Bzowski et al. 

2012; 

Mobius et al. 

20 ll). IBEX-Lo col- 


lects interstellar helium atoms over several orbits during late winter and early spring where the 
atoms are near perihelion of their trajectories ((Lee et a l.1 12012). In the earliest orbits, IBEX-Lo 
samples the “warm breeze” that consists of secondary He° atoms that are deflected beyond the 
heliopause (Kubiaket al. 2014). After that the direct interstellar He° wind dominates the count 
rate. 
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The interstellar helium wind d irec tio n was obtained from these measurements using two inde¬ 
pendent analysis methods. Mobius et al, (2012) used an analytical approach to fit the evolution of 
the helium count rates across the orbits where hejium is observed, where the outer boundary of the 
computations was at infinity. Bzowski et al. (2012) compared simulations of the helium distribu¬ 
tion obtained from the Warsaw Test Particle model with the IBEX data, and performed x 2 —fitting 
between the models and data to obtain the best helium wind parameters, where the outer boundary 
of the computatio ns w as at 150 AU, the distance at which the interaction is already underway. 
McComas et a l. (2012) resolved the small discrepancy caused by the different outer boundaries 
used in the two computations and showed that these independent analyses of the IBEX-Lo data 
yield similar results. 

Independent analysis of the 2009 and 2010 seasons of IBEX data yielded the same flow 
direction within a large acceptable range. For the IBEX flow direction during the 2009 ob- 
Bzowski et al.1 (2012) found a He flow vector of Aism,oo=79.2°, Asm,oo=-5.06°, 
and Tism, oo=6094 K. For the 2010 season, the He° flow vector was 


servmg season, 

Vism, oo=22.831 km s 1 


A ism, 00=79.2°, Asm,oo=-5.12°, Vi SM ,oo=22.710 km s \ and Ti SM ,oo= 6254 K. These solutions 
als o provi ded acceptable x 2 values for mutual fits of the four interdependent parameters (Fig. 22 


in 


Bzowsk i et al. 2012). The acceptable range for the longitude is 75.2° — 83.6°, corresponding to 


temperatures of 4400-8200 K, with each flow longitude coupled to a tightly constrained range of 
the other parameters because the uncertainties are correlated. In contrast, when the Ulysses helium 
flow longitude (Witte 2004) was used as a constraint to the IBEX longitude data, the best reduced 


X 2 fit was an order-of-magnitude larger (Fig. 17 in Bzowski et al.l 2012). Independently, and 


with the analytic method described by Lee et al. (2012), Mobius et alj (2012) determined a flow 
vector for the combined 2009-2010 seasons of Aism,oo= 79.Ot^ deg, /3ism,oo= —4.9 ± 0.2 deg, 
Vism,oo= 23.5l|o km s -1 , and Tism,oo= 6200+ 2 2oo K. The best-fitting direction of the interstellar 
helium wind using two different analysis techniques are in good agreement with each other. 


The large range of possible longitudes that is obtained when the correlated uncertainties of 
the 4-dimensional parameter are included can be constrained further if independent information 
on the LIC temperature is available. In Paper I we used the gas temperature that is obtained from 
astronomical observations of the spectra of refractory interstellar elements toward Sirius (see 


One outcome of the Bzowski et al. (2014) reanalysis of the Ulysses data set is that the “suspi¬ 
cious” coincidence remarked upon by LB 14, that the latitude of the He° flow did not vary between 
between the IBEX and Ulysses He° measurements, may not be so settled after all with the new 
reduction and modeling of the Ulysses data. Bzowski et ah (2014) found that the best-fit Ulysses 
upwind direction of the flow latitude is one degree north of the ccntctyof the IBEX range, at 
Asm,oo= 6° ± 1.0° versus /3ism,oo= 4.98° ± 1.21° for IBEX ( McComas et al. 2012 ). 
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4.1.3. IBEX instrumental response function 


In order to support their perspective that the interstellar wind direction through the heliosphere 
is invariant with time, LB 14 speculated that the instrumental response function used to analyze 
the IBEX-Lo data is incorrect, and argued that the true “dead-time” correction is 7 ms, which is 
larger than the value of the instrument response function underlying the IBEX data analysis of the 
2009-2010 observing seasons. An explanation of the incorrectness of this speculation is given in 
Appendix [0, where the details of the IBEX instrumental function are discussed. 


In addition, the ISN flow direction is driven by the variation of the peak in latitude with 
the ecliptic longitude of IBEX, which is not affected by data flow restrictions^ IBEX is a Sun- 
pointed spinner that measures vertical swaths of the sky with each spin (McComas et al. 120091) . 
The He° particle trajectories are measured over both longitude and latitude, and that are insensitive 
to the speculations of LB 14 regarding the instrumental response parameter (Appendix 14.11) . The 
best-fitting interstellar neutral (ISN) flow direction found from IBEX observations 014. 1 .21) is not 
dominated by the longitudinal height and width of the observed flux distribution. 


LB 14 have also claimed that the second independent analysis of IBEX data by Bzowski et al. 


(2012) omits “dead-time” corrections. This argument is invalid because the Bzowski et al. results 
are based on data that used an instrument response function that had been tested and validated. 


4.2. STEREO measurements of pickup ions 


Measurements of He, Ne and O pickup jons_(PUI) by the STEREO PLASTIC instrument 
between 2007 and 2011 are reported by Drews et al. (12012 . hereafter D12). Neutral interstellar 
He, Ne and O atoms penetrate to the inner heliosphere where they are ionized by photoionization, 
electron-impact ionization, or charge-exchange with the solar wind ( Bzowski et al. 2013b . § The 
charged particles that form inside of 1 AU couple to the solar wind magnetic field and are convected 
radially outwards with the solar wind where they can be observed at 1 AU. Helium and Ne survive 
the larger ionization losses associated with close passage to the Sun and become gravitationally 
focused in the downwind direction. Oxygen is mainly ionized in the upwind regions. Helium, Ne 
and O pickup ions form a region of enhanced PUI density with a crescent-like distribution in the 


9 Photoionization models of the LIC predict neutral fraction s at the heliosphere for hydrogen, helium, neon, and 
oxygen of ~ 0.78, ~ 0.61, ~ 0.20, and ~ 0.80 respectively ( Slavin & Frisch 2008 ). Up to 60%, 30%, and 4% of 
these neutral interstellar He, Ne and O atoms, respectively, survive to the Earth’s orbit before beingjonized by a mix 
of electron-impact ionization, charge-exchange with the solar wind, and photoionization (Bzowski et al. 2013). PUI 
data sample particles created between the Sun and the observation point. 
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upwind hemisphere (D12). Drews et al. determined the interstellar wind longitude independently 
from the He and Ne focusing cones, and the upwind He, Ne and O crescents, and obtained five 
independent measurements of the ecliptic longitude of the interstellar wind direction (Tabled that 
were used in Paper I to evaluate the statistical likelihood that the interstellar wind direction had 
changed over the past forty years. 

LB 14 conclude that the STEREO pickup ion data are only valid as upper limits on the in¬ 
terstellar wind longitude, and that the oxygen data should be ignored entirely. In more detail they 
argue that: (1) Pickup ion transport inside of 1 AU causes deviations of the location of the focusing 
cone and crescent of up to 5° so that values for the flow directions obtained from STEREO PUI 
can only be taken as upper limits. (2) The requirement of a smoothly varying ionization rate used 
in the STEREO analysis is not satisfied. They claim that D12 only consider photoionization and 
not considered ionization by electron impacts. (3) The inter-comparison of the data between orbit 
groupings are invalid because the data are not statistically independent. (4) The statistical simula¬ 
tions of the direction uncertainties of the STEREO analysis are incorrect. (5) Oxygen should not 
be used to infer the interstellar flow since there is also a population of secondary oxygen atoms. 


4.2.1. Transport ofPUls in the expanding solar wind 


Pickup ions are created in the turbulent environment inside of 1 AU. If the mean-free-path for 
transport of the newly-formed PUIs is large, Ay~ 1 AU, PUI ions are transported parallel to the 
magenetic field li nes so that their spatial distribution is no longer a valid marker of the ISN wind 
direction ( Chalov 2014 ). 0 However if the PUIs are formed in a non-turbulent environment, e.g. 
A||~ 0.1 AU, the parallel transport is negligible. A new MESSENGER study has shown that A ~ 
0.1 AU in the region 0.1-0.3 AU from the Sun ( Gershman et al.ll2014 ). These new results invalidate 
the main argument of LB 14 that the PUI data should not be used for determining the direction 
of the interstellar neutral wind. However, we reach this conclusion with the understanding that 
new modeling of the ionization of ISNs close the Sun is providing new insights into the stability 
of the focusing cone and crescent locations that may affect these conclusions (J. Sokol, private 
communication, 2014). 


k 


Chalov (2014J) has investigated the influence of transport effects on the helium focusing cone position for the 


unusual solar minimum conditions between 2007 and 2009, an interval that is included in the STEREO measurements. 
By numerically solving the transport equation of He + pickup ions, the possible displacement of the helium focusing 
cone was evaluated for different values of the mean free path of the atom parallel to the magnetic field, Aj|. For the 
most extreme conditions, A|| = 1 AU gave an offset angle for the observed focusing cone center of 3°, which would 
then need to be subtracted from the measured value to get the true inflowing wind direction. However for A|| =0.1 
AU of the MESSENGER data, no transport is found. 
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LB14 based their conclusions regarding the PUIs on a study ofiChalov & Fahr (120061) . who 
found that transport effects could generate a relative azimuthal displacement of PUIs to positive 
longitudes of up to five degrees, but that the transport effect was insignificant for the range of 
w = Cjon/Uw > 1-4 in the SWICS data from ACE because of instrumental properties. However 
the STEREO data used by D12 is for pickup ions with w > 1.5 so the Chalov and Fahr (2006) 
modeling does not apply. 

In addition to the new MESSENGER results, three other sets of measurements of the ISN 
flow longitude obtained from focusing cone measurements after 2007 agree with each other. The 
longitud e for the He° ISN flow found from STEREO data collected 2007-2011 is 77.4° ± 1.9° 
( Drews et al. 2012 ). It is consistent with the ACE SWICS and MESS ENGER FIPS data collected 
over the first two years of this interval. Also, Ge rshman et al. (120131 G13) analyzed ACE SWICS 
1 AU data from 2007-2009 and found a longitude 77.0° ± 1.5° that was consistent with the ISN 
flow direction, 76.0° ± 6.0°, collected by MESSENGER at 0.3 AU during three passes through the 
focusing cone during 2007-2009. Gershman et al. (2013) concluded that the similar longitudes 
found from the 1 AU SWICS data and the 0.3 AU FIPS data indicated that azimuthal transport 
effects have a minimal effect on the He° flow longitude, and that electron impact ionization is 
significant inside of 1 AU. 


However, as noted above, tr ansport effects app ear to have a minimal effect o n the direction of 
the He focusing cone. STEREO ( Drews et alil2 012). and MESSENGER and ACE (Gershman et al. 


20130 measurements of the helium focusing cone were collected conducted during similar temporal 


intervals of 2007-2009, during solar minimum, and found focusing cone directions in agreement. 
The MESSENGER data were collected at 0.3-0.7 AU while ACE and STEREO data were obtained 
near 1 AU, suggesting that transport effects are insignificant. However the larger uncertainties on 
the He° flow longitude obtained from the MESSENGER measurements, 77.0° ± 6.0°, do not allow 
possible transport effects to be ruled out conclusively. The IBEX 2009 in situ measurem ents were 


collected over a similar tim e period and agree with the PUI focusing co ne longitude (IMobius et al. 


2012 : Bzowski et al.ll2012 ). Comparisons between the ACE PUI results ( Gloeckler et al .112004 ) and 


the Ulysses in situ He° results ( Wittell2004 ). collected at similar times but different radial distances, 
identify similar helium flow directions. 


Although the possible transport of PUIs inside of 1 AU remains an important question because 
of large uncertainties on the helium focusing cone longitude in the STEREO data, and because of 
the large electron impact ionization rate that varies with time, available data suggest that it is not 
significant. 
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4.2.2. Temporal variations in PUI count rates 


The geometric properties of the gravitational focusing cone and upwind crescent vary with 
both the solar magnetic activity cycle and the solar potation. The PUIs observed at 1 AU are 
seeded by the ionization of neutrals inside of 1 AU (B zowski et al.ll2013 ). and subsequent outwards 
convection as part of the solar wind with a distinct energy distribution ( Vasvliunas & Siscoelll976 ). 


In the approach used by D12, both slowly and rapidly varying influences on the observed 
cone and crescent structure are included in the analysis, but due to the different origins of long¬ 
term and short-term modulations they are included in distinctly different ways. Slowly changing 
modulations of the pickup ion count rates include possible variations of the ionization rate induced 
by solar UV radiation or by slow decay of instrumental detection efficiency. Either would result in 
a systematic shift of the focusing cone and crescent location in ecliptic longitude. These shifts can 
go in either direction, and they depend purely on the nature of the modulation. Temporal variations 
in the interstellar wind direction were assumed to be negligible. 


Drews et ah (2012) estimated the influences of the slowly varying count rate modulation of 


interstellar pickup ions on the results and explicitly included those variations in the analysis. Com¬ 
parisons between SWICS PUI data and the solar UV emission over a ten-year period of time 
showed that the influence of slowly varying ionization rates induced by solar UV radiation is neg¬ 
ligible. Influences of a slowly decreasing detection efficiency induced by aging of the PLASTIC’S 
microchannel plates were considered explicitly in the STEREO analysis by estimating the slowly 
varying detection efficiency decay with time. Other slowly varying modulations of the observed 
pickup ion count rates, e.g. the effect of the eccentricity of the orbit of STEREO around the Sun, 
have been thoroughly discussed in D12. 


Pickup ion count rate modulation over short time scales, on the other hand, are included im¬ 
plicitly through the statistical approach used by D12. The short-term modulations mainly stem 
from solar wind density compressions, changes of the solar magnetic field, or solar wind variabil¬ 
ities such as the electron impact ionization rates (that LB 14 claim are ignored). As demonstrated 
in D12, the statistical approach is especially well-suited to deal with these ki nds of fluctuations. 


The influence of electron impact ionization at 1 AU is strongly coupled to the occurrence of 
dense solar wind streams iBzowski et al. 2013 ) and therefore fluctuates on short time scales. LB 14 
argue that the statistical analysis does not include the modulation of pickup ion count rates by 
electron impact ionization 0 This statement is misleading because the statistical analysis focuses 
specifically on using the statistics of the total data set as the basis for reliably separating out and 


1 'They comment that “Drews et al. (2012) consider only the photoionization, and omit the second source of ioniza¬ 
tion, electron impact.” 
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characterizing short-term variations such as electron impact ionizations. 

LB 14 make contradictory arguments that the analysis of the STEREO/PLASTIC data influ¬ 
ences of electron impact ionization were ignored, but also that the ISM flow direction derived 
from Gaussian fits to numerically co-added orbital counts, without futher evaluation of the data, is 
more accurate. This second method is based on pure geometry and explicitly ignores all influences 
on the foimation of pickup ions including those that are caused by variations in electron impact 
ionizations, while the D12 statistical analysis takes those influences into account. 

Further details of the analysis of the STEREO PUI data are discussed in Appendix O includ¬ 
ing the separation of temporal and spatial variations in the data set, the intercomparison between 
the four orbits acquired over the years 2007-2011, and additional statistical tests on the STEREO 
results. 


4.2.3. Effect of secondary populations 


LB 14 argue that the neutral interstellar wind (NISW) longitude that is derived from the oxygen 
PUIs in the upwind cresent is invalid because of the presence of secondary oxygen atoms. These 
atoms, however, are inconsequential for measurements of the oxygen PUIs in the upwind crescent. 

The populations of secondary interstellar neutral atoms are produced by interstellar ions 
that enter the outer heliosheath regions, and are first offset by the Lorentz force as they migrate 
through the interstellar magnetic field draped over the heliopause, and then become neutralized 
through charge-exchange or recombination so that some of them enter the heliosphere. IBEX may 


have detected the s econdary populations of He and O (M obius et al. 


2009 


Mobius et alil2012 ) but the interpretation is still not clear (llCubiak et al.ll2014h Th ese exp ected sec- 


Bzowski et al. 2012; 


ondary populations are significantly slower than the primary He and O neutrals (Izmodenov et al. 


199% . and therefore have much lower probabilities of surviving in sig nificant numbers to the in¬ 
ner heliosphere compared to the primary populations (IBzowsk i et alj 2013i) . Less than 20% of 


the interstellar oxygen is expected to be ionized (Slavin & Frisch 


20081) in the LIC, so there are 


relatively few atoms in the potential parent population for the secondary O. For a likely asymmet¬ 
ric heliosphere, secondary populations of heavy species such as oxygen will be produced away 
from the upwind-downwind axis, where their thermal speeds will be less than the plasma flow 
speed in the production region, so they are less likely to reach the 1 AU region. When the rela¬ 
tively few secondary oxygen atoms are coupled to the long travel times through the heliosphere 
for secondary populations, combined with the rapid ionization of oxyg en due to photoionization 
and charge-exchange with solar wind protons (e.g. see Bzowsk i et a l. 120131) . secondary oxygen is 
highly unlikely to alter the direction of the interstellar oxygen flow through the inner heliosphere 
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that is inferred from the O pickup ions. 


The crescent locations deduced from He + , Ne + , and 0 + pickup ion observations show a good 
agreement within the uncertainties given in D12. Because any secondary flow of neutral helium 
and neon is at best at insignificant levels in the analysis, influences by a secondary oxygen flow 
are most likely insignificant, as pointed out in D12. In addition, any effect on the 0 + crescent 
from a secondary component would push the structure to even higher longitudes than the one 
reported in D12. As [L allement et al . (120051) showed based on Lyman alpha backscatter observa¬ 
tions, the combined inflow direction of primary and secondary hydrogen is shifted towards smaller 
longitudes, which would then be expected also for the combined primary and secondary oxygen 
population. Mobius et al. (120091) indeed mentioned a structure in the neutral interstellar O flow 
that may be co nsiste nt with a secondary component shifted towards the same side as reported for 
H by Lallement et al. ( 2005 ). 


4.3. He° 584A backscattered emission 


The gravitational focusing cone was first measured in 1972 through the fluorescence of solar 
Hel 584A emission from interstellar He° close to the Sun (Weller & Meier 19741) . Since then, 
40 years of measurements of the interplanetary 584A fluorescence form a historical record of the 
interstellar wind direction ( Aiello 19781: Aiello et al. 1979 : Weller & Meier 1981 : Dalaudier et al. 


19841: IVallerga et aU2004l : lNakagawa et al.ll2008l) . The 584A data are the earliest measurements of 
the interstellar wind direction and velocity. 

LB 14 question the treatment of the He° 584A backscattering data in Paper I for three reasons: 
QJThey claim that the H c° wi nd direction used in Paper I, and obtained from the Mariner 10 data 
( Ajello 1978 : Ajello et al. 1979 ). is no t based only on t he He ° 584A emission^ (2) They contend 
that the value for the flow longitude in Dalaudier et al. (119841) that is based on Prognoz 6 measure¬ 
ments of the geometric center of the peak 584A emission from the He° focusing cone, is inaccurate 
and should not be used. (3) They argue that the_uncertainties used for the longitude determined 
from the NOZOMI 584A measurements ( Nakagawa et all 120080 are incorrect and larger than the 
value used in Paper I (±3.4°, that was taken from the original NOZOMI publication). After Paper 
I was published, an erratum for the NOZOMI measurement was published (see below), so that this 
argument of LB 14 argument was correct. It is important to note that this increase in the NOZOMI 
uncertainties leads to an statistically possible constant flow longitude of 75.0° ± 0.3° using the 
same assumptions as in Paper I, albeit at a lower likehood than a flow longitude that varies linearly 


12 They state: “We note however that the direction quoted in this work has been derived from H Lyo observations, 
namely from the downwind Lyo minimum, i.e. it applies to the neutral hydrogen and not to helium.” 
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with time (lj3] Table [2]). 


4.3.1. Mariner 10 data: 


Overlooking the discussion of the Mariner 10 data in the online supplementary material to 
Paper I, LB 14 claim that the longitude of the NISW used for the Mariner 10 data in Paper I is 
based on both He° and H° data. The resonant He° 584A emission feature was observed during two 
roll control maneuvers (RCM) o£Mariner 10 - RCM10 on 28 January, 1974 (Aiello 1978), and 
RCM3 on 8 December 1973 (Ai ello et al. 19791 . Both the He° 584A and H° Lya emissions were 
observed during RCM 10. From the simultaneous modeling of the two RCM 10 data sets, Ajello 
(1978) found the direction of the interstellar gas flow. The RCM3 observations were obtained 
while the spacecraft was located close to the axis of symmetry of the He° focusing cone so that 
the RCM3 data gave unique and independent results for the direction of the He° focusing cone. 
The RCM3 He 0 wind direction agre es w ith the wind direction in Ajello (1978) that is based on 
both He° and H° fluorescence data. 1 Ajello et alj (119791) state that a change in the location of the 
interstellar wind axis by more than 5° degraded the fits to the RCM3 584A measurements. The 
±5° uncertainties in Paper I are generous since Ajello et al. (1979) quote them for the combined 
longitude and latitude, while in Paper I we adopted the full uncertainty for the longitude alone. 


4.3.2. Prognoz 6 data 


The argument of LB 14 that the geometric wind direction from the He° 584A backscattered 
emission data from Prognoz 6 should not be used is puzzling because Bertaux is a coauthor on 
the study (Dalaudier et al. I1984h . LB 14 object to our use in Paper I of the wind direction found 
from the geometric center ofthe 584A emission pattern, even through it was quoted with uncer¬ 
tainties bv Dalaudier et alJ (119841) . The Prognoz 6 analysi s was later revised in two papers in order 


to for ce agreement with results from other experiments (IChassefiere et al.lll988l: lLallement et al. 


20041) . The meta-analysis in Paper I incorporated unbiased data published with uncertainties, so 
we assumed that both of the He° wind directions from Prognoz 6 were useful. However, in light 
of the LB 14 comments and the later revisions of the Prognoz 6 results, we provide a new fit to the 
historical data that omits both of the Prognoz 6 data points (Section |3]). 
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4.3.3. NOZOMI data 


In Paper I we used the longitude 258.7° ± 3.4° for the He° flow longitude quoted in Section 
6 of INakagawa et al.l ( 2008b . Referring to the ±3.4° uncertainties that we used for the NOZOMI 
data, LB 14 state that the NOZOMI article contains ”... no justification anywhere of such an un¬ 
certainty ...” and this is incorrect. However, LB 14 are correct with their argument that the ±3.4° 
uncertainty is too small. We now know that the ±3.4° uncertainty given in the original paper was 
a typographical error. The NOZOMI team has most kindly provided the correc ted values for the 
upwind direction of the He° flow: A = 258.78° ± 8.0°, f3 = 3.48 ± 8.0° ( Nakagawa et al. 2014 ). 
In our refit to the historical data on the He° flow direction, we have therefore adopted these cor¬ 
rected uncertainties for the He° flow longitude from the NOZOMI data. We note that the corrected 
NOZOMI longitude affects the statistical likelihood of a constant flow direction over time using 
the otherwise same data set as used in Paper I Lj3]). 


5. Conclusions 


While the LB 14 speculations about the IBEX and STEREO analyses are refuted, the differ¬ 
ences between the titles of the two papers, ”On the Decades-Long Stability of the Interstellar Wind 
through the Solar System” for Lallement & Bertaux, versus ’’Decades-Long Changes of the Inter¬ 
stellar Wind Through Our Solar System” for Frisch et al., highlights fundamental questions about 
the interaction between the heliosphere and interstellar medium that are not yet understood. 


It is shown that small-scale structure, < 330 AU, in the surrounding LIC may be present. 
Structure over these scales is indicated by the LIC collisional mean free path, possible Alfven 
waves propagating along the LIC interstellar magnetic field, and turbulence related to supersonic 
motions between the LIC and adjacent clouds that could drive shocks into the LIC. 


In response to the Frisch et al. (20131) meta-analysis of available, published data at that time, 
where it was found that the interstellar wi nd direction has more likely to have changed by ~ 6.8° 
over the past forty years than not), La llem ent & Bert aux (2014) presented a long list of arguments 
that supported their basis for concluding that the analysis was flawed. Chief among the incorrect 
speculations was that the IBEX-Lo instrumental response function is larger than the measured 
value. The suggestion that azimuthal transport of pickup ions inside of 1 AU showed that PUIs 
can only be used as upper limits on the interstellar wind direction is contradicted by results from 
the MESSENGER spacecraft on PUIs inside of of 1 AU. They also misunderstood the STEREO 
analysis procedure. The arguments of LB 14 and our rebuttal on those arguments are summarized 
in Table Q] 


The most valuable comment of LB 14 was about the uncertainties on the NOZOMI He° 
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wind direction, and it motivated our inquiry to the_NOZOMI team about those uncertainties. A 
typographical error in the original publication (INakagawa et al.1120081) has now been corrected 
( Nakagawa et alj 2014) and incorporated into our analysis Cj3]). A new statistical analysis that 
included the corrected NOZOMI uncertainty was performed on the data used in Paper I. Again, 
using the prior sets of observations, but including this correction, a temporal variation is statisti¬ 
cally indicated by the data and is still statistically highly likely. However, in contrast to the findings 
of Paper I, the larger uncertainties of the corrected NOZOMI longitude now allow a constant flow 
longitude of 75.0° ± 0.3°, although the likelihood of a constant flow is lower than the likelihood 
that the flow longitude has varied over time. The statistical fit given by header B in Table [2] should 
replace the fit that was provided in Paper 10 


The statistical fitting was repeated using improved uncertainties for the STP 72-1, Mariner 10, 
and SOLRAD1 IB results, omitting the Prognoz 6 data and using both of the 2009 and 2010 seasons 
of IBEX results. For this fit, the most likely statistical result is found from a linear fit to the temporal 
data. The variation over forty years corresponds to a longitude variation of <5 A = 5.6°±2.4°. Again, 
a constant flow longitude over time was statistically acceptable but at a lower statistical significance 
than the linear variation. We emphasize that the linear fits are not equivalent to claiming that the 
historical variation has actually occurred as a linear shift with time. 


Results from the 2012-2014 IBEX observing seasons of the flow of interstellar He° through 
the heliosphere are now being studied, however the data contains puzzling effects that suggest 
that more comp licated processes may contribute to the in teraction between interstellar He° and 
the heliosphere (McComas et al. 2014 : Leonard et ak| 2014). Those results are not included in this 
study. 


The comparison of heliospheric data on the interstellar neutral wind over long time intervals 
offers the possibility of sampling small scale structure that likely exists in the surrounding inter¬ 
stellar material. There are now over 40 years worth of historical data on the flow of interstellar He° 
through the heliosphere, and not all of it has been analyzed using contemporary knowledge of the 
He° ionization processes. Future analyses of these historical data should enrich our understanding 
of the interaction between the heliosphere and interstellar medium. 


This work was carried out as a part of the IBEX mission, funded as a part of NASAs Explorers 
Program. The researchers from SRC PAS were supported by Polish National Science Center grant 
2012-06-M-ST9-00455. The authors thank Harald Kurcharek for helpful discussions. 


13 This correction brings out the good side of the value of robust scientific exchanges. 
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A. LIC temperature from interstellar absorption lines 


Because of the cloudy nature of local ISM, 6%-33% of sp ace within 10—15 p c is filled with 
low-de nsity material that is identified primarily by cloud velocity ( Frisch et alJhoil ). Redfield & Linskv 
(120081) have identified an absorption component corresponding to the LIC velocity toward ~ 80 
stars located 2.6 to 74.5 pc away. Temperatures are reported for many of these stars, including 15 
stars that are within 22 pc (two of which are in a binary system). Temperatures are found from 
the fits to the spectral line that use an intrinsic Voight profile where the broadening is assumed to 
be due both to mass-dependent thermal broadening and a second mass-independent component. 
These temperatures are plotted in Fig. [7] against the angle between the star and the LIC downwind 
direction. The reported temperatures range between 52001 }®qo K and 120507'y^j K and show a 
tendency for the temperature to decrease toward the downwind direction. Alternatively, Redfield 
and Linsky (2008) also report a LIC temperature of 7500 ± 1300 K using a larger set of 19 stars 
out to distances of 68.8 pc. Hebrard et al. (1999) reported a LIC temperature of 8000 4 ;„j[] 0 K for 
the binary pair Sirius A and Sirius B, from the simultaneous fitting of absorption lines from both 
neutrals and ions. The range of the possible LIC temperatures plotted in Fig. [7]for the nearest stars 
indicate either that the gas temperature varies between sightlines through the LIC, or that other 
non-LIC clouds are at LIC-like velocities for some sightlines. 

In Paper I we evaluted the LIC temperature from measurements of Fe + , Mg + , and Ca + in the 
LIC component toward Sirius. The analysis was restricted to these refractory species that trace both 
ionized and neutral gasses to avoid non-overlapping distributions of neutral and ionized species 
because of the ionization gradient in the LIC (S lav in & Fri sch 120081) . A median temperature of 
58 OOI 700 K was found (see the online supplementary material of Paper I). The temperature is 
consistent with the Ulysses temperature of Witte (2004) and the prefe rred He° temperature of 
Tism,oo= 6300 ± 390 K found from IBEX data in McComas et al. ( 2012 ). This Sirius temperature 
was used as an independent constraint on the temperature uncertainties in the IBEX parameter 
tube, which then determined the flow longitude from the IBEX data that was used in the statistical 
analysis of Paper I. Temperatures derived from interstellar absorption lines are not used in the new 
fits to the data that are reported in f j3] 


B. The IBEX instrumental corrections and results for the interstellar neutral temperature 


We have determined Tism.oo from the width of the neutral flux latitude distribution in the 
analytic method Mobius et al. (2012); the value obtained through \ 2 -minimization by Bzowski et 
al. (2012) is also dominated by the latitudinal variation of the fluxes. The peak height of these 
distributions, and thus the width, could possibly be affected by a suppression of the data flow of 
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individual events through the data system that depends on particle count rate. To account for such a 
potential effect, at that time still unknown, Mobius et al. (2012) estimated an upper limit in terms of 
a maximum event processing time of r < 5 ms. They modeled the suspected effect as an instrument 
dead-time and compared observations of the same ISN flow at IBEX-Lo energy steps 2 and 4, the 
latter with a rate lower by a factor of 5 and thus reduced interface traffic. Lallement & Bertaux 
(2014) have argued that an instrumental dead-time r ~ 7 ms would allow a reconciliation of the 


angular distributions observed with Ulysses and IBEX, although iMobi u s et al. (2012) showed that 
r < 5 ms is an upper limit and McComas et al. (2012) reported that the effect on the temperature 
is much smaller based on follow-up testing and simulation, as discussed in more detail below. 

The telemetry limitation for high-flux neutral atom sources, such as the ISN flow and the 
Earth’s magnetosphere, was anticipated, and energy-angle histograms of H and O with 6° reso¬ 
lution are accumulated onboard IBEX to normalize the reported event rates. Because the width 
of the IBEX-Lo field-of-view is already 7° FWHM and to eliminate any potentially unanticipated 
effec ts on the data flow by the limited telemetry, the angular resolution used by Bzowski et al. and 
Mobius et al.1 (2012) was restricted to 6° 


The hardware interface between IBEX-Lo and the data handling system has a known trans¬ 
fer time of 1.354 ms for each individual event, which increases for the first 1° portion of each 
6 ° bin to 2.188 ms because housekeeping data and count rates are transmitted after each event. 
Both transmission times are much shorter than the previously estimated maximum dead-time. In 
addition, the interface has a double buffer, which reduces the resulting suppression even further. 
The response of the interface to a large variety of event and rate combinations was tested with a 
hardware simulator and computed with Monte Carlo simulations. These tests and reviews of the 
data system hardware and software rule out response times longer than the nominal values quoted 
above. 


It was found that the majority of the event traffic across the interface, of typical rates of 300- 
400 s^, orig inate s fro m an unanticipated substantial electron flux through the sensor, as mentioned 
by Mobius et al. (2012). Electrons are identified based on their short time-of-flight and thus are 
completely eliminated from ENA signals. However, these events add an almost isotropic (or slowly 
varying) base rate that burdens the interface. Figure [6] shows simulation results for the combination 
of such a base rate and an angular ISN flow distribution, with a peak rate of 75 events/s (the max¬ 
imum ISN event rate observed) as a function of the base event rate. Shown are the total absolute 
reduction factor for the base and peak event rate combination (open symbols), the reduction factor 
for the peak relative to the already reduced base rate (solid black), the resulting increase of the 
width of the angular distribution (solid red), and the increase of the deduced temperature (solid 
blue). While the total event flow is reduced for the observed base rates of 300-400 events/s by 
up to 20%, the relative reduction of the peak rate over the base is only 2.5-3%. Only the latter 
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reduction is relevant for any distortion of the angular distribution, and consequently the width is 
increased only by ~ 1%, which translates into an apparent increase in the deduced temperature of 
only 2%. Interestingly, these effects on the relative height of the peak distribution and the deduced 
temperature top out for a base rate of ~ 700 events/s at < 4% for the peak rate and < 2.5% for 
the temperature. Any systematic instrumental effects on the temperature determination from IBEX 
observations larger than 2% and on the peak fluxes as a function of longitude larger than 2.5% can 
be excluded for the observations analyzed thus far. 

It should be noted that the type of events causing this uncertainty are the unanticipated elec¬ 
trons that are of very little use in the IBEX ENA analysis, and they have been excluded by a new 
instrument mode that started mid-2012. This operational modification now eliminates these reduc¬ 
tion effects almost entirely. A direct comparison of angular ISN flow distributions in both modes 
has verified in flight the total event rate reduction shown in Figure [5] and demonstrated that no 
significant broadening of the angular distribution is seen when the base event rate is present. 


C. Separating the temporal and spatial variations in the STEREO data 


The analysis strategy was explained in detail in Drews et al. (2012). D12 disentangle the 
temporal and spatial variations in the STEREO data through the use of equation 3 in their paper, 
where the pickup ion count rates are separated into the terms representing the time-independent 
production model for the pickup ions, D r (</>) for the azimuthal angle 0, and the time-dependent 
modulation parameter M(f), which describes the fluctuations of pickup ion count rates induced by 
solar wind variabilities on short time scales. Under the assumption that the production mechanism 
is the same for each orbit, an approximation of the modulation parameter, G(t), can be determined 
from the variability of the pickup ion count rates with respect to the average of all of the orbits, as 
is described through equation 4 in their paper. The possible effects of short-term variabilities of 
high-speed solar wind streams or electron impact ionization were overcome by fitting the product 
of the time-dependent modulation parameter, G(t), and a model for the angular distribution of 
interstellar pickup ions A) to the angular distribution of He + , 0 + , and Ne + observed by 

STEREO/PLASTIC (equation 6 in their paper). Not including G(t) into the fitting procedure, i.e., 
G(t) = constant, would mean disregarding the temporal fluctuations of pickup ion count rates 
induced by solar wind variabilities such as electron impact ionization, and is therefore comparable 
to a standard fitting technique, i.e., fitting a Gaussian distribution to the focusing cone and the 
pickup ion crescent. 


LB 14 have also objected to the method of comparing data between two separate orbits using 
different possible permutations of the four orbits, arguing that the different orbit combinations “are 
not fully independent”. We disagree with this viewpoint. Due to the random nature of the pickup 
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ions fluctuations in each orbit, the intrinsic comparisons between pairs of different orbits produces 
independent information. The dependent variable that arises from these different permutations of 
orbit comparisons is the invariant component of the pickup ion distribution that does not change 
with time and depends only on the inflow direction of the interstellar helium atoms. The indepen¬ 
dent variable is the estimate of the relative modulation parameter G(t), which was deduced from 
equation 4 in their paper and changes for each permutation of orbit combinations. Because G(t) 
is a key parameter to recover the inflow direction from pickup ion data, the intercomparison of 
the individual orbits is required to determine the uncertainty of the modulation parameter, A G(t), 
which in turn is necessary to determine a scientifically justified estimation of the uncertainty of the 
inflow direction of interstellar matter, A A. 


Further simulations of the STEREO data are presented in Drews (2013). This extension of 
the statistical analysis compares two different approaches, the numerical co-adding of the count 
rates that ignores the time variability of the STEREO data, versus simulations of the properties 
of each orbit compared to the overall data sample. These comparisons were performed for several 
different sets of conditions and are presented in Drews (2013), chapter 4.3. The statistical approach 
was based on 1500 randomly generated data sets using four different kind of modulations. Drews 
compared the results of the simulations to those of a standard fitting technique. Drews (12013l Fig. 
4.4) showed that the standard fitting technique provides a less precise fit to the STEREO data, 
producing an uncertainty for the focusing cone and crescent position that is ~2 times larger than 
the uncertainty found using the full statistical analysis. 
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Table 1. Comments on Lallement & Bcrtaux ( 2014 ) criticisms of data and usage in Paper I 



Argument 

Comment 

Section 

(Appendix) 

1 

Criticisms regarding IBEX analysis: 
Ulysses and IBEX velocity-longitude 

True for hot models using Maxwellian 

14.1.114.1.21® 

2 

combinations give same perihelion 
IBEX-Lo instrumental response 

distributions, but irrelevant to results 

Incorrect; the IRC, dead-time correction, 

irnidm 

3 

correction (IRC) is wrong 

Correct IRC makes IBEX 

and data throughput are well understood 
Incorrect; the IBEX results rely on 

ECU 

4 

and Ulysses longitudes agree 

Warsaw modeling based on 

variations of the peak as a function of 
latitude and observer longitude 

Incorrect; Warsaw models compared to 

EU 

5 

wrong IRC corrections 

Criticisms regarding STEREO analysis: 
Short-term electron-impact ionization 

data based on a correct and tested IRC 

Incorrect; short-term Eli is distinguished 

0021 

6 

(Eli) not included in data analysis 

Orbit groupings are not 

from long-term ionizations in methodology 
Incorrect; each orbit pairing selects 

0(0 

7 

statistically independent 

Simulations of uncertainties of 

invariant distributions from comparisons 
of data acquired over separate time-spans 
Uncertainty simulations are 

0 © 

8 

STEREO data are incorrect 

Secondary oxygen invalidates usage 

appropriate for variable PUI data 

Incorrect; secondary oxygen population 

E23] 

9 

of oxygen PUIs 

Longitudes from PUI focusing cone 

minor compared to primary population, and 
if included push longitudes into the 
opposite direction than suggested 

Maybe, but transport effects appear to be 

EZH 

10 

and upwind crescent are upper limits insignificant, based on comparisons between 
due to transport effects inside 1 AU several different data sets 

Criticisms regarding use of Ultraviolet Data: 

Mariner 10 longitude is based on Incorrect: Mariner 10 papers quote 

ELD 

11 

both H° and He° 

He° flow direction from Prognoz 6 

longitude of He° wind separately 

The Prognoz 6 He° direction from geometric 

021 

12 

geometric analysis not valid 

NOZOMI uncertainties of ±3.4° are 

arguments has uncertainties and was included 
since it was not repudiated in the literature 
True; the results of the 2014 NOZOMI 

E33] 


too small 

erratum are used in this paper 



A Numbers in last column show section or appendix where the topic is discussed. 
































Table 2. Fits to historical interstellar wind data (March 3, 2015) 


Num. Basis 

of fit 

Fit* 

A(deg.) 

Reduced 

x 2 

p-value 

A. Paper I : A : 




1 Linear 

70.6(±1.6) + 0.17(±0.06) * t 1970 

0.97 

0.49 

2 Constant A 

75.1(±1.3) 

1.71 

0.031 

B. Corrected NOZOMI uncertainties ; B 



3 Linear 

71.3(±1.6) + 0.14(±0.06) * fi97o 

0.7 C 

0.82 

4 Constant A 

75.0(±0.3) 

1.1 

0.35 

C. Two IBEX data points, modified 584A uncertainties ; D 



5 Linear 

71.4(±1.6) + 0.14(±0.06) * *1970 

0.8 

0.3 

6 Constant A 

75.0 ±0.3 

1.3 

0.2 

7 Parabolic A fit 

74.1(±1.4) - 0.11(±0.11) * *1970 ± 5.0(±2.1)10 -3 * f? 970 

0.6 

0.14 


* A is the ecliptic longitude of the interstellar wind and £1970 is the elapsed time since 1970. A p-value larger 
than 0.05 passes the p-test and therefore is not unlikely. 

A This fit is the original linear fit in Paper I that used the NOZOMI value from Nakagawa et al. (2008). 

B These fits used the historical wind data in Paper I, except that the uncertainties for the NOZOMI value are 
corrected to incorporate the erratum of Nakagawa et al. (2014). This linear fit should replace the fit listed in Paper 
I. See Section l4~3l 

c This fit is accepted and is more likely than fit number 4. The small value of \ 2 suggests that the uncertainties 
may be overestimated by ~ 20 %. 


Nakagawa et al. (2008 


D These fits incorporate uncertainties for the wind longitude for the STP2-1, Mariner 10, and SOLRAD1 IB data 
that include only the longitude uncertainties (see Section [431) . The geometric Prognoz 6 direction is omitted from 
the fit according to the discussion in Section 14.31 The He° wind directions from the individual IBEX 2009 and 
2010 observing seasons are included independently (Sections 14.1.21 [31 and see Bzowski et al.ll2012 ). 
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Fig. 1.— Left: The fit using the correct NOZOMI uncertainties is shown by the dashed line, and 
the dotted line shows the original (Paper I) fit with the uncorrected NOZOMI uncertainties. The 
shaded region shows the uncertainties for the new fit. The uncertainties on the IBEX data point 
(red dot in upper right) are constrained by the IBEX parameter range and the LIC temperature 
toward Sirius based on Fe + , Mg + , and Ca + data (Paper I, OSM-S4). Right: Same fit as on the 
left, except that the Prognoz 6 data are omitted, the uncertainties for the STP72-1, Mariner 10, 
and SOLRAD11B data are revised to include only the longitude uncertainty, and seasons 2009 
and 2010 of the IBEX data are listed as independent data points. The IBEX longitude in the right 
hand figure does not depend on any constraints that are based on temperatures obtained from LIC 
absorption lines. These fits are explained further in Section [3} 
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Fig. 2.— Hyperbolic trajectory geometry in the plane of the ISN bulk flow. The location of the 
perihelion (or flow tangential to the Earth’s orbit) is degenerate in all combinations V ISMi 00 (Ai S m,oo) 
that satisfy equation (1). The solid and the dashed blue line indicate two sample trajectories from 
infinity to their perihelion at 9 = 0° that satisfy the relation, but arrive from different inflow 
longitudes Aism.oo as indicated by the solid and dashed black lines. Aobs is the longitude of the 
observations, 9^ is true anomaly of the inflow trajectory (blue line) from infinity to perihelion at 
r E = 1 AU. The dashed blue line shows the trajectory of an atom with a different #oo-Vism,oo 
combination that satisfies equation (1). 
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Fig. 3.— Composite of angular distribution measurement as a local Mach Cone, whose width 
(shown on the left) = f t h(iAU)/^iSM(iAU) is defined by the ratio of the local thermal and bulk 
velocities of the flow. The angular width in latitude a ^ is obtained from the cut through the maxi¬ 
mum of the neutral helium rate distribution in ecliptic longitude (along the white arc) shown in the 
color coded map (center panel). The latitudinal rate distribution of that cut (right panel) is shown 
along with a Gaussian fit and statistical uncertainties. 
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Fig. 4.— Schematic view of the interstellar gas flow in its trajectory plane (blue), which is inclined 
to the ecliptic (yellow) at the latitude inflow angle /3ism,<x>- V’Peak is the latitude of the flow in 
the solar rest frame, 6^ is the perihelion angle with respect to the interstellar flow at infinity, and 
Asm,oo is the latitude of the interstellar flow at infinity. The x,y,z (and primed) coordinates define 
the cartesian coordinate systems. 
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12 14 16 18 20 

Orbit 


Fig. 5.— Peak location in latitude (top) and width (bottom) of ISN flow distributions for the IBEX 
orbits during the 2009 ISN flow season, simulated for the actual observation periods with two ISN 
flow parameter sets along the center of the IBEX tube in parameter space. Thcj inflow longitude 
of Witte (2004) is shown in red and the one of Bzowski et al. (2012) and Mobius et al. ( 2012 ) 
in blue. The latitudinal width of the distribution, a, is also shown in green for the temperature 
found by Witte (2004). These simulations are based on the measured properties of the IBEX-Lo 
instrument. A temperature of 8180 K would be required for the Witte (2004) flow vector to match 
the latitudinal width measured by IBEX. 
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Fig. 6.— Simulated effects on the observed flux distribution of the ISN flow for the maximum peak 
rate of 75 events/s. Absolute total rate reduction (open symbols), reduction of the ISN flow peak 
rate relative to the base rate reduction (solid black), increase of the angular width a in latitude (solid 
red), and resulting apparent temperature increase (solid blue) as a function of different continuous 
base event rates. 
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Fig. 7.— The temperatures of the LIC toward stars within 22 pc that sample the LIC, according to 
Table 1 of iRedfield & Linskvl (120081, dots), are plotted against the angle between the star and the 
ISN downwind direction {i = 184°, b = —12° in galactic coordinates). The triangle shows the 
LIC temperature from Paper I, which was determined from the IBEX parameter range combined 
with Fe + , Mg + , and Ca + absorption line data in the LIC component toward Sirius. The wide 
green shaded region shows the_ uncertainty on the temperature range found from the reanalysis of 
the Ulysses data by Bzowski et alj (2014). The blue and dark-green bars show, respectively, the 
temperatures required by IBEX and Ulysses data if the true interstellar flow upwind direction is 
the best-fit value found from the Ulysses data, or 255.3° according to B zowsk i et alJ (2014). The 
difference between these two temperatures (853 K) shows that IBEX and Ulysses data are not 
in agreement. The dashed line shows a second order fit through the interstellar temperature data 
(dots), and is meant to guide the eye and illustrate the systematic decrease in reported cloud tem¬ 
perature as the angle to the downwind direction decreases. If this apparent temperature decrease is 
spurious, it could represent the misidentification of LIC components in sidewind directions. 



























